Mach cone shock waves at RHIC 

JORG RUPPERT^ 

Department of Physics, McGill University, 3600 University Street, Montreal, QC, 

H3A 2T8 Canada 

AND 

Thorsten Renk 

Department of Physics, PO box 35 FIN-40014, University of Jyvaskyla, Finland, 
and Helsinki Institute of Physics, PO Box 64, FIN-00014, University of Helsinki, 

Finland. 



Energy and momentum lost by hard jets propagating through hot and 
dense nuclear matter have to be redistributed in the medium. It has been 
conjectured that collective sound modes are excited. Those lead to Mach 
cone nuclear shock waves in the nuclear medium that are shown to account 
for three and four particle angular correlation structures of hadrons with a 
(semi-)hard trigger hadron in heavy-ion collisions at RHIC. 

1. Introduction 

The quenching of QCD jets created in relativistic nuclear collisions has 
been proposed as an important indicator of the creation of a quark-gluon 
plasma [IJ. It is extensively studied theoretically and experimentally at 
RHIC. The main emphasis in most theoretical studies (before 2005) has 
been solely on the description of the radiative energy loss which the lead- 
ing parton suffers while traversing the nuclear medium due to the emis- 
sion of partonic secondaries. In this paper we focus on another aspect of 
the in-medium jet physics, namely the question if the jets traversing the 
medium can transfer energy and momentum to collective modes in the nu- 
clear medium that might be able to account for the emergence of peculiar 
signals in the particle correlation measurements. Recently, measurements 
of two- and three-particle correlations involving one hard trigger particle 
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have shown a surprising sphtting of the away side peak for all centralities 
but peripheral collisions, qualitatively very different from a broadened away 
side peak observed in p-p or d-Au collisions [2]. Interpretations in terms 
of colorless [3l, '3] and colored [5j sound modes have been suggested for an 
explanation of this phenomenon. For an overview also discussing alternative 
mechanisms, see [6] and references therein. While the microscopic excita- 
tion mechanism of the colorless modes is under investigation [7] , a detailed 
theoretical study of the experimentally observable correlations is already 
now possible if this mechanism is assumed to effectively excite the mode. In 
the following, we discuss how such Shockwaves lead to observable correlation 
signals in the dynamical environment of a heavy- ion collision. 



2. The model 

A Monte Carlo simulation of the hard back to back process in the 
medium is performed. There are four main stages in the modeling: 1) 
the primary hard pQCD process, 2) the description of the soft medium, 3) 
the energy loss from hard to soft degrees of freedom, 4) the simulation of 
the Shockwave propagation and its modification of the soft medium until de- 
coupling. For brevity we only present a sketch of the model here, a detailed 
description of the model can be found in [8] . The soft medium is described 
by a parametrized evolution model [9] which gives a good description of the 
bulk matter transverse momentum spectra and HBT correlation radii. The 
energy loss for a given parton path inside this medium is described prob- 
abilistically by P{AE, E)pa,th, which is the probability for a hard parton 
of energy E to lose energy AE while traversing the medium in the ASW 
formalism [10] . Since we are interested in the energy deposition on aver- 
age in a given volume we focus on the distribution of the average energy 
< AE >= P{A{E))AEdAE along the paths to infer dE/dx, see Fig. 1 
in [8] . We assume that a fraction f of the energy lost to the medium excites 
a Shockwave characterized by a dispersion relation E = CsP where Cg is the 
speed of sound inferred from the equation of state by = dp/de. The 
dispersion relation determines the initial angle of propagation of the shock 
front as (j) = arccosc^. The time dependent energy transfer to the mode per 
unit time is given by fdE/dr. Each piece of the front is propagated with 
the local speed of sound through the medium. Once a wavefront has evolved 
until it satisfied freeze-out condition T = Tp it cannot propagate further 
leading to an additional boost at freeze-out. The Cooper-Frye formalism is 
employed to calculate the hadronic distribution accordingly. On the near- 
side the trigger condition is realized. The important role of longitudinal 
and transverse expansion as well as flow effects on the observed correlation 
signal are discussed in [8]. 
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3. Results 

In order to show how the excitation of a collective mode can account 
for the observed correlation signal in a quantiative manner, we show in 
Fig. [1] the correlation signal and its sensitivity to different descriptions 
of the soft background. We chose an energy independent / = 0.75 which 
accounts for the observed signal. We show results for two different expansion 
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Fig. 1. Correlation signal on the away side for different soft background medium 
density distributions and evolutions (see text, not acceptance-averaged). 

patterns (Bjorken and non-Bjorken evolution) and two different transverse 
densities (box and nuclear profile Ta)- Due to the position of the freeze- 
out hypersurface at large radii for the box density, the shockwave gets on 
average longer exposure to transverse flow in this scenario which leads to 
a less pronounced maximum since transverse flow can erase a peak if flow 
and shockwave are not aligned. The angle is larger for longitudinal Bjorken 
expansion which is due to the fact that initial cooling is rapid and therefore 
the averag temperature quickly approaches the phase transition temperature 
Tc where Cg is small. In this specific calculation we suppressed a full rapidity 
averaging, leading to a somewhat larger angle than expected if averaging 
were performed. In Fig. [2] we show the resulting 3- particle correlation for 
non-Bjorken evolution with initial box density and nuclear density profile, 
for details of the calculation see [H]- The region close to the trigger is not 
studied. Our results fit nicely with the fact that transverse flow can erase 
in some configurations one or even both of the wings of the cone. In this 
case, only the diagonal region is populated whereas if both wings appear 
the off diagonal maxima are populated. In this sense 2- and especially 3- 
particle correlation measurements do not only reveal information about the 
smallness of Cs in the evolution, indicating the existence of a (cross-over) 
phase transition to a quark gluon plasma, but also can eventually be used 
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to further constrain our understanding of the evolution of the soft-matter 
at RHIC. This work was supported by the Natural Science and Engineering 
Research Counsil of Canada and the Academy of Finland. 
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